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The effect of surface topography on the nonlinear 
dynamics of Rossby waves 

By S.I. Abarzhi, 0. Desjardins, H. Pitsch 


1. Motivation and objectives 

Boussinesq convection in rotating systems attracts a sustained attention of the fluid 
dynamics community, because it has intricate non-linear dynamics (Cross & Hohenberg 
1993) and plays an important role in geophysical and astrophysical applications, such as 
the motion of the liquid outer core of Earth, the Red Spot in Jupiter, the giant cells in 
the Sun etc. (Alridge et al. 1990). A fundamental distinction between the real geo- and 
astrophysical problems and the idealized laboratory studies is that natural systems are 
inhomogeneous (Alridge et al. 1990). Heterogeneities modulate the flow and influence 
significantly the dynamics of convective patterns (Alridge et al. 1990; Hide 1971). The 
effect of modulations on pattern formation and transition to turbulence in Boussinesq 
convection is far from being completely understood (Cross & Hohenberg 1993; Aranson 
& Kramer 2002). 

It is generally accepted that in the liquid outer core of the Earth the transport of the 
angular momentum and internal heat occurs via thermal Rossby waves (Zhang et al. 
2001; Kuang & Bloxham 1999). These waves been visualized in laboratory experiments 
in rotating liquid-filled spheres and concentric spherical shells (Zhang et al. 2001; Kuang 
& Bloxham 1999). The basic dynamical features of Rossby waves have been reproduced 
in a cylindrical annulus, a system much simpler than the spherical ones (Busse & Or 
1986; Or & Busse 1987). For convection in a cylindrical annulus, the fluid motion is 
two-dimensional, and gravity is replaced by a centrifugal force, (Busse & Or 1986; Or & 
Busse 1987). Hide (1971) has suggested that the momentum and heat transport in the 
core might be influenced significantly by so-called bumps, which are heterogeneities on 
the mantle-core boundary. To model the effect of surface topography on the transport 
of momentum and energy in the liquid outer core of the Earth, Bell & Sowarcl (1996), 
Herrmann & Busse (1998) and Westerburg & Busse (2001) have studied the nonlinear 
dynamics of thermal Rossby waves in a cylindrical annulus with azimuthally modulated 
height. 

The models of Bell & Sowarcl (1996), Herrmann & Busse (1998), and Westerburg & 
Busse (2001) considered different regimes of the flow, and reported analytical and numer- 
ical solutions for the system of conservation laws. Bell & Soward (1996) assumed periodic 
boundary conditions in the radial direction and neglected the viscous friction in the fluid 
interior. They found two types of convective patterns depending on the modulation am- 
plitude: Rossby waves for small modulation amplitude and traveling waves with long 
azimuthal length scale for large modulation amplitude. Herrmann & Busse (1998) and 
Westerburg & Busse (2001) have accounted for the viscous friction in the fluid interior, 
neglected the friction produced by the Eckman layer, assumed no-slip boundary condi- 
tions in the radial direction, and found Rossby waves with quasi-periodic dependecies 
in space and in time. Both rigorous studies, the work by Bell & Soward (1996), and by 
Herrmann & Busse (1998) and Westerburg & Busse (2001) compliment each other and 
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identify several important features of the modulated convection. A need still remains for 
a model, which quantitatively and qualitatively describes the effect of modulations on 
Rossby waves in a wide range of the forcing parameters and convection intensity, and 
which can compare this inhomogeneous system to other pattern-forming systems (Cross 
& Hohenberg 1993; Aranson & Kramer 2002). Here we suggest a model to study the 
effect of surface topography on the nonlinear dynamics of thermal Rossby waves. 


2. Pattern-forming system with modulated forcing 

For rotating cylindrical annulus with homogeneous boundary conditions, the nonlinear 
dynamics of thermal Rossby waves are described by a complex Ginzburg-Landau equation 
in terms of amplitude A , whose magnitude determines the convection intensity, and whose 
phase describes changes in the position and direction of the convective rolls. We refer 
the reader to Busse & Or (1986) and Or & Busse (1987) for details. For heterogeneous 
systems, a rigorous derivation of the amplitude equation from the conservation laws has 
not been accomplished yet (Cross & Hohenberg 1993; Aranson & Kramer 2002). To 
model the nonlinear dynamics of slightly inhomogeneous systems, one usually applies 
the Ginzburg-Landau equation with forcing. As an example, one may consider pattern 
formation in the presence of a mode associated with a conserved physical quantity, such as 
the Rayleigh-Benard convection driven by surface tension, which has been considered by 
Siggia & Zippelius (1981), Tribelsky & Velarde (1996), Matthews & Cox (2000). Other 
examples are convection in liquid crystals driven by a stochastic electric field (Meyer 
et al. 1987; Wu & Andereck 1990; Roder et al. 1997), and the study of the effect of 
time-periodic forcing on an oscillatory system (Elphick et al. 1987; Coullet et al. 1990; 
Hemming & Kapral 2000). 

Convection in a rotating cylindrical annulus with inhomogeneous boundaries differs 
from the foregoing pattern- forming systems. The surface modulation is not a slowly 
relaxing Goldstone mode as in the models with a conserved physical quantity (Siggia 
& Zippelius 1981; Tribelsky & Velarde 1996; Matthews & Cox 2000). The modulation 
does not result in an Ising-Bloch bifurcation as in the model with time-periodic forcing 
(Elphick et al. 1987; Coullet et al. 1990; Hemming & Kapral 2000), and its amplitude 
does not change with time, in contrast to the system with stochastic forcing (Meyer et al. 
1987; Wu & Andereck 1990; Roder et al. 1997). To study the effect of surface topography 
on the dynamics of thermal Rossby waves, the Ginzburg-Landau equation with spatially 
periodic and time-independent forcing is a proper model. The equation has the form: 

^ + AoZl + D\A\ 2 A + (Se ikY + S*e~ ikY ) A . (2.1) 

Here T and Y are the slow time and coordinate, A is the complex amplitude, the values 
of f3 and D are complex numbers with (3 = (3 r + ify, D = D r + iD.i, the subscripts r and 
i mark the real and imaginary parts, Ao ~ ( R — R cr ) is real with R cr being the critical 
Raleigh number (Cross & Hohenberg 1993; Busse & Or 1986; Or & Busse 1987; Aranson 
& Kramer 2002). The forcing amplitude and the absolute value of the wave- vector are 6 
and k respectively, where the asterisk marks the complex conjugate. The system (2.1) is 
slightly inhomogeneous and <5|fc <C 1. 

Without forcing, (5 = 0, the values of /?, D , A 0 in the equation (2.1) can be derived 
from the incompressible Navier-Stokes and energy equations for R ~ R cr (Busse & Or 
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1986; Or & Busse 1987). In agreement with the stability criterion the real parts of /3 and 
D obey the relations (3 r > 0 and D r < 0 (Cross & Hohenberg 1993; Aranson & Kramer 
2002). The dynamics of the convective rolls are descirbed by the traveling waves with 
A = A = Aoe zqv+llT , where 

Al = - (A ° , 7 = -ft rf + D.A,. (2.2) 

This solution appears for Aq > /3 r q 2 , and depicts the pattern formation of thermal Rossby 
waves without forcing. 

If the value of A 0 is finite and the modulation amplitude is small, |d| <C | A 0 1 , the 
explicit dependence on the Y coordinate in equation (2.1) can be eliminated. We present 
the amplitude A in (2.1) in the form A = A+A, where A = Aoe lqy+l ' lT from (2.2), and the 
term A = A1+A2+... appears due to the forcing with A n = (a n e~ mkY + b n e mky ) e lqy + * 7T , 
n = 1, 2, .... In this way we obtain 

A = e iqY+i ' rT [A 0 + ai e~ ikY + b ie ikY + a 2 e~ 2ikY + b 2 e 2ikY + ...] . (2.3) 


Substituting the expression (2.3) in (2.1) and expanding the equation (2.1) in terms of 
small |<5/Ao|, we can derive a system of the coupled equations of the Ginzburg-Landau 
type: 


dA a 

dT 


[/o + 7 ?Aq A 2 2D ( a i®i A b-^b* + fli^i)] Aq + «i<5* + bi S , 


(2.4) 


da 1 

W 


fi Q\ -\~dUi T [D (aia^ T 2&x^i) a\ T 2 DAq (a^a2 A a 2 bi -t- bib2')]AAQS-\-a 2 5 (2.5) 


dbi 

dT 


da\-\- f Y bi~\~[D (2aia* + did*) b\ + 22? Ag (aia^ + a\b 2 + bib 2 )]-{-Ao5* -\-b 2 S, (2.6) 


do 

— f 2 a 2 A db\ A [2?Ao ( 2a±bi + af)] + aid + 03 d* , (2-7) 

db 

—— = da\ + f ~2 b 2 + [2?Ag (2a*6i + b\ )] + did* + 63 d ... (2-8) 


where /o = Ao — (3q 2 — ij, f± = Ao — (3 (q ± nk ) 2 — *7 + 22? Ag, d = DAq , and the 
Y -dependence of the values A 0 , a n , b n has been neglected for the sake of simplicity. 

The nonlinear dynamics in the system (2. 4-2. 8) depend on the forcing parameters and 
the convection intensity, with the modulation wave-vector k being the key factor for 
pattern formation. Below we briefly describe several limiting cases. 

In the linear approximation in | d / Ao | <C 1, the solution (2.2) satisfies equation (2.4), 
while the equations (2. 5, 2. 6) are reduced to: 


da\ 

W 


db 1 
dT 


fi o-i A db\ + Agd, 


da\ + f+ bi + A 0 d . 


(2.9) 
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The linear system (2.9) has steady solutions with 


A 

O 

* 

1 

AoS(ff-d*) 

(fr (ft)* - dd*) ’ 

{fr(ftr-M-) 


(2.10) 


The solutions (2.2,2.3,2.10) are stable for Re[a ] < 0, where the Liapunov exponents a 
obey the equation 


(cr + 2/ 0 ) (a - f 1 ) (a - (ft)*) - dd* 


= 0 . 


( 2 . 11 ) 


For the large-scale modulations with k/q <C 1, pattern formation is quite complicated. 
In the interval /3 r q 2 < A 0 < F r , the pure solution with A 0 from (2.2) and cq = bi = 0 
occurs, where F r = F r {D , /3 ) and F r > (3 r . For A 0 > F r q 2 , this solution becomes unstable, 
and the modulated pattern (2.3,2.10) appears. In the limit of {k/q) —> 0, the values of 
| o-i | , 1 6i | increase as |<5/Ao|(<7/fc) 2 , and result in a growth of the higher-order terms in 
the system (2. 5-2. 8). A nonlinear coupling among the amplitudes a n and b n yields then 
a solution with \a n /A 0 \, |&„/A 0 | |<5/A 0 |, yet keeping |a n /A 0 | <C 1. In this case the 

system (2. 5-2. 8) is reduced to 


da 

grj-, = fi d\ + dbi + [ D (did± + 2b\bi ) d± + 2DAq (a\d 2 + d2b\ ^ 1 ^ 2 )] 1 (2.12) 

db 

= d a i + fi^i + [D (2aia\ + bib±) bi + 2 DAq ( aid \ + ai &2 + 6 *^ 2 )] ? (2.13) 

dd 

— fi a 2 + d &2 + \DA$ (2aibl + a 2 )] , (2-14) 

db- 

= dat + f 2 b 2 + [DAq (2a*6i + 6 2 )] . (2.15) 

From (2.12-2.15) we obtain 

a n = A n e inGT , b n = B n e inGT , n= 1,2,3, (2.16) 


with A 2 , B 2 dependent on Ai = Ai{A 0 ,k,q) and B\ = Bi{A 0 ,k,q). Using the solv- 
ability conditions for real A\, B i, in (2.12-2.15), we derive the dispersion relation G = 
G{A 0 ,k,q). It is remarkable that G ~ 7 for small q and for {k/q) <C 1. The solutions 
and their stability analysis are quite cumbersome and not presented here. 

For the small-scale forcing with k/q 1, the modulated state (2.3,2.10) is stable for 
all A 0 > 0. The nonlinear dynamics in this case are governed by terms proportional to 8 
in (2. 5-2. 8 ), and the nonlinear solutions are steady. With a 2 = b 2 = 0, equations (2. 4, 2. 5) 
are reduced to 


dai 

tff 


Aq8 T D (qu* T 2b\bf) q , 


(2.17) 
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Bb 

=A 0 S* + D(2a 1 al + b 1 bt)b 1 , (2.18) 

and for small but finite |<5/Ao| in the system (2.17,2.18) 

bi ~ A 0 \5/Aq\ 2 ^ . (2-19) 

We conclude that the large-scale forcing eventually results in the occurrence of Rossby 
waves with quasi-periodic time-dependence. The convective pattern (2.3,2.16) consists 
of three traveling waves with different group velocities. The dynamics of this pattern 
depends on the convection intensity and the forcing wave- vector, yet it is independent of 
the forcing amplitude. In contrast, the small-scale forcing produces spatially modulated 
convective patterns with dynamics governed by the forcing amplitude. 


3. Numerical results 

To obtain numerical solutions of the complex amplitude A in (2.1), we non-dimensionalize 
this equation first using the length of the domain L, the reference time scale A)) 1 , and 
the initial magnitude of A 0 . This leads to the following equation: 

f) A' Pi 2 A' 

fyffi = gyrl + ^ + D'\A'\ 2 A' + S' A' (3.1) 


The initial solution of the problem has the form A = Aoe lqY , and the forcing is 
defined as S = 6e lkY + 5*e~ xk ^ . The non-dimensional quantities are given by A! = A/Aq , 
Y' = Y/L , T 1 = TA 0 , P' = P/A 0 L 2 , D' = DAq/A 0 and S’ = 5/ A 0 , and with the 
wavenumbers k and q being k' = kL and q' = qL. 

We solve equation (3.1) on a domain of length L = 10 with periodic boundary con- 
ditions. We replace the complex equation by a set of two equations for the real and 
imaginary parts of A. This set of equations is discretized using a second order central 
difference scheme in space. A fully implicit solver is used in time. 

For these computations we choose A 0 = 100 and A 0 = 1. The values of P and D 
are taken to match experiments of Herrmann & Busse (1998); Westerburg & Busse 
(2001). The three remaining parameters are the forcing amplitude S and the ratio be- 
tween the two wavenumbers k and q. Several cases have been computed with S /Aq = 
0.1, 0.2, 0.4, 0.6, 0.8, 1, and ( k/q ) = 1/10, 1/3, 1, 3, 10. The space is discretized with a 300 
points mesh, and the time step is chosen to be AT = 0.05. 

Figure 1 presents the real part of A' in (3.1) for the case with no forcing. The pattern 
has the form of the traveling wave, and its dynamics agrees qualitatively and quantita- 
tively with the solution (2.2). 

Figures 2, 3, and 4 show the real part of the amplitude, and the spatial and temporal 
Fourier transforms in the case of large-scale forcing with k/q = 1/3 and 5/Aq = 0.2. 
Three traveling waves with different group velocities, described in (2.19), are clearly 
visible on the Figure 2. The spatio-temporal characteristics of the pattern are in good 
agreement with the results in (2.16). Figure 3 shows that the asymptotic dynamics of 
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Figure 2. Real part of the amplitude in the case of large-scale modulation 


the pattern are highly nonlinear, and there exist three major spatial modes associated 
with the wavevectors q, q — k and q + k. Temporal Fourier transforms of the pattern on 
presented Figure 4 clearly show the existence of a low frequency, associated with small 
(7 — G) <C 7 in (2.16). The other modes on Figure 4 correspond to 7 in the equation 
(2.2) and (7 + G) ~ 27 in the solution (2.16). We conclude that the large-scale forcing 
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Wavenumber 

Figure 3. Spatial Fourier transform of the amplitude in the case of large-scale modulation 


x 10 4 



Frequency 


Figure 4. Temporal Fourier transform of the amplitude in the case of large-scale modulation 


influences the Rossby waves significantly, and forms a pattern modulated in space and 
quasi-periodically in time. 

Figures 5, 6, and 7 present the real part of the amplitude, and the spatial and temporal 
Fourier transforms in the case of small-scale forcing with k/q = 10 and 6/ A 0 = 0.1. The 
pattern formation in this case differs significantly from that of the large-scale modula- 
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Figure 6. Spatial Fourier transform of the amplitude in the case of small-scale modulation 


tions. Figures 5, 6, and 7 clearly show the formation of spatially modulated Rossby waves, 
in agreement with the solution (2.19). Figure 6 demonstrates that the spatial dynamics 
of the pattern are described by the modes associated with the wavevectors q, q — k and 
q + k. Temporal Fourier transforms of the pattern shown on Figure 7 indicate that there 
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Figure 7. Temporal Fourier transform of the amplitude in the case of small-scale modulation 


is only one frequency in the system, which is associated with 7 in (2.2). We conclude 
that the small-scale forcing modulates Rossby waves in space, (2.19). 


4. Discussion 

In the present work, we studied the effect of surface topography on the nonlinear 
dynamics of thermal Rossby waves. The foregoing results show that the interaction of 
Rossby waves with spatially modulated forcing gives rise to a quasi-periodic behavior 
in space as well in time. The large-scale forcing results in the appearance of interacting 
traveling waves, whose dynamics depend strongly on the forcing wavevector and weakly 
on the convection intensity, and are insensitive to the forcing amplitude. In contrast, the 
small-scale forcing generates spatially modulated patterns, but does not influence their 
temporal dependence. In this case, the pattern dynamics are governed by the forcing 
amplitude and the convection intensity. 

Our model captures the main properties of the linear and weakly nonlinear solutions 
obtained in Bell & Soward (1996) and Westerburg & Busse (2001). On the other hand, 
the model describes the influence of modulations on the dynamics of Rossby waves in 
a wide range of the forcing parameters and convection intensity, and predicts the new 
properties of the nonlinear dynamics. Our results are in a good qualitative agreement 
with the observations by Westerburg & Busse (2001). 
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